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Abstract. Precise measurement of neutrino oscillations, and hence the determination of 
their masses demands a quantitative understanding of neutrino-nucleus interactions. To 
this aim, two-body meson-exchange currents have to be accounted for along within real¬ 
istic models of nuclear dynamics. We summarize our progresses towards the construction 
of a consistent framework, based on quantum Monte Carlo methods and on the spectral 
function approach, that can be exploited to accurately describe neutrino interactions with 
atomic nuclei over the broad kinematical region covered by neutrino experiments. 


1 Introduction 

The description of neutrino interactions with nuclei, besides being interestir^ in its own right, provides 
an essential input for neutrino-oscillation experiments, such as DUNE yj], which plan to measure 
neutrino oscillation parameters and the neutrino mass hierarchy and determine the charge-conjugation 
and parity (CP) violating phase. Such experiments make use of nuclear targets as detectors that, while 
allowing for a substantial increase in the event rate, demand a quantitative understanding of neutrino- 
nucleus interactions. 

The MiniBooNE collaboration has reported a measurement of the charged-current quasielastic 
(CCQE) neutrino-carbon inclusive double differential cross section, exhibiting a large excess with 
respect to the predictions of the relativistic Eermi gas model ||2]- It has been suggested that the dis¬ 
crepancy might be due to the occurrence of events with two particle-two hole final states Hi that 
are not taken into account by the relativistic Eermi gas. Within a realistic model of nuclear dynam¬ 
ics, their occurrence arises naturally owing to two-body meson-exchange currents and/or correlations 
induced by the nuclear interaction. 

Because neutrino beams are always produced as secondary decay products, their energy is not 
sharply defined, but broadly distributed. Therefore, the observed cross section for a given energy 
and angle of the outgoing lepton includes contributions from energy- and momentum-transfer regions 
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Figure 1. Schematic representation of the inclusive electron-nucleus cross section at beam energy around 1 GeV, 
as a function of energy loss, from Ref. (^. 


where different mechanisms are at play. They can be best identified in electron-scattering experiments, 
in which the energy of the incoming electron is precisely known. The typical behavior of the double 
differential inclusive cross section for the process 

£? -b A —> ^ ^ (1) 

with a beam ener^ around 1 GeV where only the outgoing electron is detected, is shown in Figure[Tl 
taken from Ref. 151]. The target nucleus in its ground state and the undetected hadronic final state 
are denoted by A and X, respectively. At small energy transfer the structure of the low-lying energy 
spectrum and collective effects are important. In the quasielastic peak region, the cross section is 
dominated by scattering off individual nucleons although a significant contribution from nucleon pairs 
is also present. The A-resonance region is characterized by one or more pions in the final state, while 
at large energy transfer the deep inelastic scattering contribution is the largest. 

The ground state of the nucleus does not depend upon the momentum transfer of the process and 
can be accurately described within nonrelativistic many-body theory (NMBT). On the other hand, the 
final state does depend on the momentum transfer and the nonrelativistic approximation can be safely 
applied only in the region of low and moderate momentum transfer, corresponding to jqj < 500 MeV. 
The treatment of the region of high momentum transfer requires a theoretical approach in which the 
accurate description of the nuclear ground state provided by NMBT is combined with a relativistically 
consistent description of both the final state and the nuclear current. 

Setting up a framework in which two-body currents and nuclear correlations induced by realistic 
Hamiltonians are consistently accounted for is the main topic of this proceeding. 


2 Lepton-Nucleus cross section 


The double differential inclusive cross-section of the electromagnetic (EM) and neutral-current (NC) 
lepton-nucleus transition {+A —> {' +Xm proportional to the product of the leptonic and the hadronic 
tensors. 
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( 2 ) 
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The leptonic tensor is completely determined by lepton kinematics, whereas the hadronic tensor, 

W,y = 2 <017/|X)<X|7y|0> 6^^\po + q-px), (3) 

X 

containing all the information on strong interaction dynamics, describes the nuclear response of the 
target to the transition currents 7^ from the initial state |0), with four momenta po, to the final state 
IX), with four momenta px- 

The nuclear current consists of one- and two-nucleon contributions; the latter arise from processes 
in which the interaction with the beam particle involves a meson exchanged between the target nucle¬ 
ons and are usually referred to as meson-exchange currents (MEC) 

7- = + j; y = Xi tr + Z Ju + Z t', ^ (4) 

i j>i k>j>i 

The relation between the MEC and the nuclear hamiltonian will be discussed in the following section. 

3 Moderate momentum transfer regime 

In the moderate momentum transfer regime, both the initial and final states appearing in Eq. Q are 
eigenstates of the nonrelativistic nuclear Hamiltonian H 

//|0)-£o|0) , H\X)^Ex\X). (5) 

The nuclear Hamiltonian describing the dynamics of the point like nonrelativistic nucleons is given 
by 


/4 2 7 ^ ^4 

"“ZSi+Z'w* z <« 

/=1 7>i=l k>j>i=l 

In the above equation, p, is the momentum of the i-th nucleon, while vij and Vijk are the two- and 
three-nucleon potentials, respectively. 

Eor nuclei as large as the ground state wave function can be obtained via the Green’s Eunction 
Monte Carlo (GEMC) approach ||a,|3l- It has to be remarked that, within the limits of applicability of 
NMBT, the GEMC is truly an ah initio approach, allowing exact calculations of a number of nuclear 
properties. 

In the moderate momentum-transfer regime, the nuclear cross section can be rewritten in terms 
of the response functions R^viq, a>), obtained from Eq.Q replacing the components of the current 
operator with their expressions obtained in the relativistic limit. 

A key feature of the description of neutrino-nucleus interactions at low and moderate momentum 
transfer is the possibility of employing a set of electroweak charge and current operators consistent 
with the Hamiltonian of Eq.®. The nuclear electromagnetic current, 7em = (7gj„, Jem)Arivially related 
to the vector component of the weak current, is constrained by H through the continuity equation |@] 

V ■ Jem + i[H, 7/j] = 0 . (7) 

Since the two- and three-nucleon potentials Vij and Vijk do not commute with the charge operator 7/,, 
the above equation implies that 7em involves two- and three-nucleon contributions, as shown in Eq.®. 

The one-body electroweak operator is obtained from a nonrelativistic expansion of the covari¬ 
ant single-nucleon currents, while the two-body charge and current operators that we employ in our 
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Figure 2. Euclidean electromagnetic transverse response function of aX. q = 570 MeV (adapted from Ref. 
11811 1. The theoretical predictions obtained including only one-hody and both one- and two-body transition op¬ 
erators are represented by open circles and solid circles, respectively. The shaded band refers to the Euclidean 
response extracted from Ref. ( 13 . 


calculations are derived within the conventional meson-exchange formalism Nonrelativistic 

MEC have been used in analyses of a variety of electromagnetic moments and electroweak transitions 
of s- and p-shell nuclei at low and intermediate values of energy and momentum transfers lEMl, 
improving on the description of the experimental data with respect to the one-body approximation. 

We have used the GFMC method to calculate the imaginary-time response function-the so-called 
Euclidean response function-to the electromagnetic and neutral weak currents of ^He and '^C. The 
Euclidean response function is defined as the Laplace transform of the response, 

J f>CO 

dtoe^™Rafi(q,aj), ( 8 ) 

^th 


where aith is the inelastic threshold and the Ca/x are q'-dependent normalization factors. In Ra/siq, w) the 
tu-dependence enters via the energy-conserving d-function and the dependence on the four-momentum 
transfer - q^ - aP' of the electroweak form factors of the nucleon and At-to-A transition in the 
currents. Once the latter dependence has been removed, as described in Ref. jlSll . the Euclidean 
response can be expressed as a ground-state expectation value, 


Ea/x(q,T) ^ <0|(9l(q)g-t"-^°i^O/;(q)|0) 
Ca/xiq) (0|e-(«-®o)nO) 


In the above equation H is the nuclear Hamiltonian of Eq. (|6]l (here, the AV18H-IL7 model), t is the 
imaginary-time, and Eq is a trial energy to control the normalization. The calculation of the above 
matrix element is carried out with GFMC methods 123 similar to those used in projecting out the 
exact ground state of H from a trial state. 

It has to be noted that evaluating the matrix element of Eq. da for a nucleus as large as 
is computationally nontrivial, as it requires computer capabilities and resources available on only 
forefront computers. 

In figure |2] taken from Ref. 11811 . the electromagnetic transverse Euclidean response function 
of '^C, is compared to the one obtained from the analysis of the world data carried out by Jourdan 
in represented by the shaded band. Two-body MEC contributions substantially increase the one- 
body electromagnetic transverse response function. This enhancement is effective over the whole 
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Figure 3. Electromagnetic transverse response functions of ''He at ^ = 500 MeV. Experimental data are from 

Ref. m. 


imaginary-time region we have considered, with the implication that excess transverse strength is 
generated by two-body currents not only at energies larger than the one corresponding to the quasi¬ 
elastic peak, but also in the quasi-elastic and threshold regions. The full predictions for the transverse 
Euclidean response functions is in excellent agreement with the experimental data. 

The inversion of a Laplace transform subject to statistical Monte Carlo errors, needed to retrieve 
the energy dependence of the responses, is long known to involve severe difficulties. However, there 
are techniques developed in condensed matter theory and other contexts, that seem to have success¬ 
fully overcome the inherent ill-posed nature of the problem. One of these is known as the maximum 
entropy technique |21, we have recently used it to perform stable inversions of the "^He electro¬ 


magnetic Euclidean response. 

The "^He electromagnetic transverse response function (at q = 500 MeV), obtained from inversion 
of the corresponding Euclidean response is shown in figure (31 The inversions are, to a very large 
degree, insensitive to the choice of default model response IlSll needed by the maximum entropy 
method. Results obtained with one-body only (dashed line) and (one-Hwo)-body (solid line) currents 
are compared with an analysis of the experimental world data 12311 (empty circles). There is excellent 
agreement between the full theory and experiment. Two-body currents significantly enhance the trans¬ 
verse response function, not only in the dip region, but also in the quasi-elastic peak and threshold 
regions, providing the missing strength that is needed to reproduce the experimental results. 


4 Relativistic regime 

The dynamical model discussed in section [3 can be also employed to describe the nuclear response 
in the kinematical region of large momentum transfer, in which relativistic effects play a major role. 
In this regime, the final state includes at least one particle carrying a large momentum ~ q, and fully 
relativistic expressions need to be retained for the transition currents. 

The impulse approximation (lA) scheme and the spectral function formalism allow one to cir¬ 
cumvent the difficulties associated with the relativistic treatment of the nuclear final state and current 
operator, while at the same time preserving essential features (such as correlations) inherent to the re¬ 
alistic description of nuclear dynamics described in the previous section. The lA scheme is based on 
the fact that in the kinematical region of large momentum transfer, in which |q|“' «c d, with d being 
the average nucleon-nucleon separation distance in the target, the nuclear scattering can be reason- 
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ably assumed to reduce to the incoherent sum of elementary scattering processes involving individual 
nucleons. 

The ground state of the target nucleus, which does not depend on q, is still an eigenstate of the 
nuclear Hamiltonian of Eq. |6l Within the lA the final state factorizes 

|X> = |p)(8>|«a-i,P«>- (10) 


In the above equation, the state |p) is a plane wave describing a noninteracting nucleon, while 
Iua-i, Pn) is an eigenstate of the (A - 1) nuclear Hamiltonian of Eq. (|6]l describing the (A - l)-system 
recoiling with momentum p„. Using the above equation, it is possible to rewrite the nuclear transition 
matrix element of the one-body current between free nucleon states - which can be computed exactly, 
retaining the fully relativistic expressions of the currents - and the nuclear amplitude involving the 
target ground state and the state of the spectator (A - l)-system. By using Eq. ([Tol l in Eq. Q, the 
nuclear cross section can be written in the form J^l 

da-iA^J^ / d^kdEPi(k,E)da-i, (11) 

where do-j is the cross-section describing scattering on the individual i-th nucleon, the momentum and 
removal energy of which are distributed according to the spectral function P,(k, E). Highly accurate 
theoretical calculations of the spectral function can be carried out for uniform nuclear matter, exploit¬ 
ing the simplifications arising from translation invariance The results of these calculations have 
been combined with the information obtained from coincidence (e, e'p) experiments at moderate en¬ 
ergy and momentum transfer, to obtain spectral functions of a variety of nuclei within the local density 
approximation (EDA) 13 ■ 

Note that, because of the strong spatial-spin-isospin correlations that are present among nucleons 
in the nucleus, the recoiling (A - l)-nucleon system is not necessarily left in a bound, one hole, state. 
In fact, two hole-one particle states, in which one of the spectator nucleons is excited to the continuum, 
typically contribute 15 - 20% of the spectral function normalization, the corresponding strength being 
located at large momenta and energies. 

Neglecting the contributions of final states involving more than two nucleons in the continuum, 
the cross section can be written as 


dcr = dcr.pih + dcr2p2h . 


( 12 ) 


Carrying out calculations of nuclear amplitudes combining fully relativistic MEC currents and a 
description of nuclear dynamics taking into account short range correlations requires a generalization 
of the factorization ansatz. The authors of Refs. 1 3 3 have extended the factorisation ansatz of Eq. 
(ITOt to allow for a consistent treatment of the amplitudes involving one- and two-nucleon currents. 
The resulting expression is 

|2f) = |PP'> ® |»lA-2, Pm) , (13) 


where the states |pp') and |OTA- 2 ,Pm) describe two non-interacting nucleons of momenta p and p' 
and the (A - 2)-particle residual system, respectively. The two-particle two-hole contribution to the 
hadronic tensor obtained from the extended factorisation ansatz is a sum of three terms. 


W, 


-2p2h 




2p2h 


/IV 11 


W2p2h ^ ^2p2h ^ 
fiv 12 fiv22 ■ 


(14) 


where j* and W ^^22 involve the squared amplitudes of the matrix elements of one- and two- 

nucleon currents, while [2 describes the interference between the amplitudes involving one- and 
two- body currents. 
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Figure 4. Electromagnetic transverse response of *^C at momentum transfer {q| = 570 MeV, taken from Ref. m. 
The solid line represents the results of the full calculation, whereas the dashed line has been obtained including 
only amplitudes involving the one-body current. The contributions arising from the two-nucleon current are 
illustrated by the dot-dashed and dotted lines, corresponding to the pure two-body current transition probability 
and the interference term, respectively. The experimental data are taken from Ref. (T^. 


Figure m taken from Ref. OTT . shows the transverse electromagnetic response of at |q| = 
570 MeV computed using the carbon spectral function of Ref. ll25n and approximating the two-hole 
spectral function of carbon with that of uniform nuclear matter at density corresponding to Fermi 
momentum kp =221 MeV. The fully relativistic expression of the MFC described in Refs. 1^, 29] 
was used, with the same form factors and A-width. The solid line represents the results of the full 
calculation, whereas the dashed line has been obtained including only the amplitudes involving the 
one-body current. The contributions arising from the MFC are illustrated by the dash-dot and dotted 
lines, corresponding to the pure two-body current transition probability and the interference term, 
respectively. The latter turns out to be sizable, its contribution being comparable to the total two-body 
current response for w < 350 MeV. It has to be noted that these results still need to be improved, as 
they do not include the corrections taking into account the effects of final state interactions. The data 
resulting from the analysis of in are also included for comparison. 

A comparison between the results of figure 0] and the GFMC results of figure [3 shows distinctive 
discrepancies in both magnitude and energy dependence of the two-body current contributions. While 
part of the disagreement is likely to originate from using fully-relativistic and nonrelativistic MFC, as 
well as from the non relativistic nature of the GFMC calculations, the large interference contribution 
in the region of the quasi elastic peak observed in figure [3 may well arise from interference between 
amplitudes involving the one- and two-body currents and Iplh final states. 
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